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Reactive dynamics of inertial particles in nonhyperbolic chaotic flows 
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Anomalous kinetics of infective (e.g., autocatalytic) reactions in open, nonhyperbolic chaotic flows are im- 
portant for many applications in biological, chemical, and environmental sciences. We present a scaling theory 
for the singular enhancement of the production caused by the universal, underlying fractal patterns. The key 
dynamical invariant quantities are the effective fractal dimension and effective escape rate, which are primarily 
determined by the hyperbolic components of the underlying dynamical invariant sets. The theory is general 
as it includes all previously studied hyperbolic reactive dynamics as a special case. We introduce a class of 
dissipative embedding maps for numerical verification. 
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Many chemical and biological processes in fluids are char- 
acterized by a filamental distribution of active particles along 
fractal invariant sets of the advection chaotic dynamics. These 
fractal structures act as dynamical catalysts for the reaction, 
which is relevant for a variety of environmental processes in 
open flows, such as ozone depletion in the atmosphere |1|| and 
population dynamics of plankton in the oceans 12']. 

The study of active processes in open chaotic flows has 
attracted a great deal of interest from the dynamical system 
community fsf I?]. Most of the studies have been performed 
for time-dependent two-dimensional incompressible flows, in 
the limit of weak diffusion. The flow is nonturbulent, but 
the particle dynamics is considered to be chaotic (Lagrangian 
chaos) and the active particles to interact with one another 
without modifying the flow. The advection dynamics of such 
particles can be cast in the context of chaotic scattering, where 
incoming tracers spend some time in a mixing (scattering) re- 
gion before being scattered along the unstable manifold of the 
chaotic saddle. As a result, the products of the reaction con- 
centrate along a fattened-up copy of the unstable manifold, 
giving rise to the observed fractal patterns. 

Although filamental patterns have been observed in nature, 
a clear relation between the observed value of the fractal di- 
mension and the underlying advection dynamics has been 
lacking. For example, in the "flow past a cylinder" system 
previously considered ||3|, the dimension of the unstable man- 
ifold is known to be 2 but the relevant dimension governing 
infective and collisional reactions is about 1.6. This lack of 
relation, while not reducing the importance of the previous 
phenomenological characterization of filamental distributions 
of active particles, has led to some skepticism about the merit 
of the dynamical system approach to the problem. In general, 
the advection dynamics can be characterized as either hyper- 
bolic or nonhyperbolic. In hyperbolic chaotic scattering, all 



the periodic orbits are unstable and there are no Kolmogorov- 
Arnold-Moser (KAM) tori in the phase space, while the non- 
hyperbolic counterpart is frequently characterized as having 
both chaotic and marginally stable periodic orbits. Fundamen- 
tal assumptions in such works are: (I) the active particles are 
massless point-like tracers and (2) the advection dynamics of 
these particles is hyperbolic iQl- However, in realistic situa- 
tions, the Lagrangian dynamics is typically nonhyperbolic and 
the active particles have finite size and inertia. Indeed, fully 
hyperbolic systems are quite rare and represent very idealized 
situations as the advection dynamics of tracers in fluids is usu- 
ally constrained to have a nonhyperbolic character because of 
no-slip boundary conditions at the surface of obstacles. Ob- 
stacles are at the same time the origin of Lagrangian chaos and 
the origin of nonhyperbolicity. Even away from obstacles and 
boundaries, chaotic motions of tracers typically coexist with 
regular motions. In addition, the individual active particles are 
often too large to be regarded as noninertial, as is the case for 
many species of zooplankton in the sea. Therefore, nonhyper- 
bolic and inertial effects are prevalent in nature and expected 
to play an important role in most environmental processes. A 
question of physical importance is then: what happens to the 
reactive dynamics when assumptions (1) and (2) are dropped? 

In this article, we present a scaling theory for the reactive 
dynamics of inertial particles in nonhyperbolic chaotic flows. 
The key concepts in our framework are the effective fractal 
dimension and effective escape rate, which are respectively 
defined as 
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where N{e) is the number of e-squares needed to cover the 
relevant fractal set, and R{n) is the fraction of particles that 
takes more than n ~ n{e) steps to escape from the mixing 
region (see below). As a representative application of these 
concepts, we show, for autocatalytic reactions of the form A + 
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B 2B, that the area covered by _B-particles in the steady 
state obeys the following scaling law: 
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where /ie// = (^e// + k)t, t is the time interval between 
successive reactions (time lag), tr is the reaction range, and 
k is the contraction rate due to dissipation. For nonhyper- 
bolic flows in two dimensions, D^f / < 2 and n^f / > are 
non-trivial functions of the scale e, which can often be re- 
garded as constants over a wide interval, even though D — 
Van^^o Def fie) = 2 and k, = \\m.^->Q k^j f{e) = 0. We 
find, surprisingly, that D^f / and Kef / are significantly differ- 
ent from D and n, respectively, not only for noninertial but 
also for inertial particles, even though the advection dynamics 
of the latter is hyperbolic, meaning that scars of the nonhy- 
perbolic conservative dynamics are observable in the hyper- 
bolic dynamics of slightly dissipative systems (k ^ 1). The 
previous relations for noninertial particles in hyperboUc fluids 
appear as a particular case of our results. 

The nature of the chaotic scattering arising in the context of 
particle advection in incompressible fluids may change fun- 
damentally as the mass and size of the particles are increased 
from zero. Physically, this happens because of the detachment 
of the particle motion from the local fluid motion. For spher- 
ical particles of finite size, the particle velocity v = dx/dt 
is typically different from the (time dependent) fluid velocity 
u = u(x, t) and, in first order, is governed by the equation [6] 
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The parameters are a = 3pf/{pf + 2pp) and a — ^a/St, 
where pf and pp are the densities of fluid and particle, re- 
spectively, and St is the Stokes number, which goes to zero 
in the limit of point-like particles |7|. For neutrally buoy- 
ant particles, the mass ratio parameter is a = 1, while for 
aerosols and bubbles, we have a < 1 and a > 1, respec- 
tively. The inertia parameter a determines the contraction rate 
or dissipation in the phase space (x, v), which for incom- 
pressible flows can be shown to be —2a. In the limit a ^ oo, 
the dynamics is projected on a surface defined by v = u, 
which corresponds to the advection dynamics of point par- 
ticles. The configuration-space projection of the particle mo- 
tion is strongly influenced by a 1 8 ] . For small inertia (large a), 
V ■ V w a-i(a - 1)V ■ ((v • V)u) = a'^a ~ l)(s2 - u:'^), 
where s and uj are proportional to the strain rate and vortic- 
ity of the fluid f^, respectively. The behaviors of bubbles 
and aerosols are then qualitatively different. For instance, 
along a closed orbit, aerosols are pushed outward, while bub- 
bles are pushed inward. We first consider bubbles, whose 
configuration-space dynamics is dissipative when the vorticity 
overcomes the strain rate. 

Dynamically, the inertial effects are effectively those due to 
dissipation, so that the transition to finite inertia is equivalent 
to a transition from open Hamiltonian to dissipative dynamics. 
It has been recently shown LI 0.1 that, while hyperbolic dynam- 



ics is robust, nonhyperbolic chaotic scattering typically un- 
dergoes a metamorphosis in the presence of arbitrarily small 
amount of dissipation. For nonhyperbolic scattering in open 
Hamiltonian systems, particles can spend a long time in the 
neighborhood of KAM tori, resulting in an algebraic decay 
for the survival probability of particles in the scattering re- 
gion. As a consequence, the fractal dimension of the invariant 
manifolds is the phase-space dimension |11]. This should be 
contrasted with the hyperbolic case, whose decay is exponen- 
tial and fractal dimension is typically smaller. The dissipation, 
however, may convert marginally stable periodic orbits of the 
KAM islands into attractors. The survival probability then 
becomes exponential, the dimension of the chaotic saddle be- 
comes fractional, and the overall dynamics of the scattering 
process becomes hyperbolic. 

To understand the meanings of D^f / and n^f f for fractal 
sets arising in the transition from Hamiltonian nonhyperbolic 
to weakly dissipative chaotic scattering, we consider a Cantor 
set, which is constructed in the interval [0, 1] according to the 
rule that in the nth time step, a fraction A„ = -f / {(3 + n) + S 
is removed from the middle of each one of the N = 2"^^ re- 
maining subintervals, where 7, and 6 are constants. The 
conservative case corresponds to (5 = 0, which is charac- 
terized by an algebraic decay with n of the total length re- 
maining, given by R{n) n^^ for n ;» /3, and by a unity 
fractal dimension for the invariant set, D = 1. The re- 
moved fraction A„ decreases at each time step and, as a re- 
sult, a systematic change of scales is induced, resulting in a 
non-self-similar invariant set that becomes denser as we go 
to smaller scales. The relevant consequence is that the box- 
counting dimension converges slowly to 1, leading to a scale- 
dependent effective fractal dimension, -De// ~ 1 — 
for small e. Similarly, the effective escape rate behaves as 
Keff ~ 7/n ~ 7ln2/lne^^, where n = n{e) is defined as 
the number of iterations needed to make the length of each 
remaining subinterval smaller than e. 

The limiting dynamics changes drastically and acquires 
properties of hyperbolic dynamics when a small amount of 
dissipation is allowed, which is modeled by < (5 7//?. 
In particular, the total length of the remaining intervals de- 
cays exponentially, R{n) ~ (1 — 5)" for n ^ — fi, 
and the dimension of the invariant set becomes smaller than 
one, namely, D = In 2/ In [2/(1 — S)]. At finite time, how- 
ever, the transition from the conservative to the dissipative 
case is much smoother. For Ine^^ ^ (3, the effective frac- 
tal dimension and the effective escape rate are Deff{e) ~ 
ln2/ln[2/(l - 5)]~y/lne-^ and Keff{e) « ln(l-(5)-i + 
7' hi [2/(1 — S)] / In e^^ , respectively, where 7' = 7/(1 — 6). 
The key feature is that unrealistically small scales are required 
to resolve the limiting values of the fractal dimension and 
the escape rate, rendering them physically irrelevant. For in- 
stance, to obtain -De// > 0.95, scales e < lO^'^" may be 
required. Thus, the physically important characteristics of the 
fractal set are the effective dimension and escape rate. 

We now present a physical theory for the scaling law Q, 
valid for autocatalytic reactions in two-dimensional time- 
periodic flows. Consider the area covered by B-particles in 
the open part of the flow (i.e., region where particles eventu- 
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ally escape to infinity) and, to be specific, that the time lag t is 
integer multiple of the flow's period. After a sufficiently long 
time from the onset of the reaction, the reagent B is distributed 
along stripes of approximately uniform width, mimicking the 
unstable manifold. The average width e of the stripes changes 
aperiodically over time until the steady state is reached, when 
it undertakes the periodicity r of the reaction. We assume 
that the reaction is sufficiently close to the steady state so that 
Deffie) and Keff{f-) can be considered constant over time. 
This condition is not so restrictive because for many systems 
De/ / and n^f f are essentially constant over several decades 
(see below). Therefore, for scales larger than e, the area cov- 
ered by _B-particles can be regarded as a fractal characterized 
by dimension Deff{e) and escape rate Ke//(e). 

Let e("-i)(T) and e(")(0) denote the average widths of the 
stripes right before and right after the nth reaction 1 3], respec- 
tively. Between successive reactions, the stripes shrink due to 
escape and dissipation as follows: e(")(r) = e(")(0)e~'''=^^ 
where h^jf = {neff + ^)/(2 — Df,ff) plays the role of an 
effective (contracting) Lyapunov exponent, while k accounts 
for the nonconservative contribution. When the reaction oc- 
curs, the widening due to the reaction is proportional to the 
reaction range: e'"+-'^'(0) — e*^"-'(T) cx a. The area covered by 
_B-particles right before the (n + l)th reaction satisfies 



A 



These relations can then be com- 



bined to yield the following recursive relation for the area: 
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where 

{Keff + 'i)T and c is a constant geometric factor 



From the condition A 



= A'^g^ follows our main scal- 
ing (|3} for the area Ab in the steady state 1 12]. This scaling 
holds for both noninertial and inertial particles, regardless of 
whether the flow is hyperbolic or nonhyperbolic. The hyper- 
bolic case with inertial particles is studied in Ref. fl3ll . The 
scaling (|3} represents a further step toward generality since it 
is also valid for nonhyperbolic flows. 

To make possible a numerical verification of the scaling law 
(|3}, it is necessary at present to use discrete-time maps. To 
construct a class of maps that captures all essential features of 
continuous-time chaotic flows, we note: (1) the fluid dynam- 
ics, determined by dx/dt = u, is embedded in the particle's 
advection equation and is recovered in the limit a — > cxo; (2) 
the phase-space contraction is determined by a (irrespective 
of a); (3) for small inertia, the configuration-space contrac- 
tion is proportional to (a — 1). For an area preserving map 
x„_|_i = M(x„), representing the dynamics of a time-periodic 
incompressible fluid, a possible choice for the corresponding 
embedding map representing the inertial particle dynamics is 
x„+2 - M(x„+i) = e^"(Q!X„+i - M(x„)), where the fac- 
tors involving a and a are naturally imposed by the particle 
dynamics I14ll . This can be written as 



x„+i = M(x„) + <5„, 

Sn+i = e"°(Q:x„+i - M(x„)), 



(5) 
(6) 



where x and 5 can be interpreted as the configuration-space 
coordinates and the detachment from the fluid velocity, re- 
spectively, so that (x, 5) represents the phase-space coordi- 



nates. This class of embedding maps can be a paradigm to 
address many problems in inertial advection dynamics as it 
captures the essential properties of Eq. In particular, it is 
uniformly dissipative, with phase-space contraction rate equal 
to e^^"; the noninertial dynamics x„+i = M(x„) is recov- 
ered in the limit a ^ oo; and the configuration-space contrac- 
tion rate is proportional to e^"'{a — 1) for e"°(a — 1) ^ 1, 
in agreement with the distinct behavior expected for aerosols 
and bubbles. Therefore, for finite a, a rich higher dimensional 
dynamics with a-dependent x-space projection is expected. 
Next we consider such a dynamics for both a > 1 and a < 1. 

To simulate the flow, we consider a two-dimensional area- 
preserving map that has a pronounced nonhyperbolic char- 
acter 111]: {x,y) — > [X{x — w"^ / A) , {y + w^)], where 
w = X + y/A and A > 1 is the bifurcation parameter The 
dynamics is nonhyperbolic for A < 6.5. For A 4, for exam- 
ple, there is a major KAM island in the xy-space, as shown 
in Fig.^a). Also, from Fig.^a), one can see tangencies be- 
tween the stable and unstable manifolds in the neighborhood 
of the KAM island, which is a signature of nonhyperbolic- 
ity. It is well established that, within the nonhyperbolic re- 
gion, the dimension of the invariant manifolds is 13 = 2 and 
the escape rate is k = ITTil . When this map is em- 
bedded in Eqs. ( I5l6t . for a > 1, the xy-projection of the re- 
sulting 4-dimensional map is dissipative in the mixing region 
(KAM islands and neighborhood). In this regime, the dissi- 
pation stabilizes marginally stable periodic orbits in the KAM 
islands of the conservative map, converting the KAM islands 
and neighborhood into the corresponding basin of attraction of 
the newly created attractors, as shown in Fig.^^b). The basin 
itself extends around the mixing region mimicking the stable 
manifold of the conservative dynamics. As a result, the tan- 
gencies between the invariant manifolds apparently disappear, 
suggesting that the advection dynamics of bubble particles is 
hyperbolic. For a < 1, on the other hand, the configuration- 
space projection expands in the mixing region and almost all 
the orbits eventually escape to infinity. However, for small in- 
ertia and a close to 1, particles in the regions corresponding to 
KAM islands of the conservative dynamics and neighborhood 
are almost trapped in the sense that the time it takes to escape 
is much larger in these regions than outside them. These re- 
gions are neglected in our analysis of the open part of the flow, 
as shown in Fig. [He), because filamental structures cannot be 
resolved inside them. 

Numerical simulation of the autocatalytic reaction is per- 
formed by dividing the mixing region with a grid where the 
size of the cells represents the reaction range a. Particles are 
placed in the center of the cells. When a reaction takes place 
in a cell occupied by a B-particle, all the cells adjacent to 
it are infected with B-particles 1.15.1 . We assume that A is 
the background material and that the reaction takes place si- 
multaneously for all the particles at time intervals r fl^. If 
we start with a small seed of B-particles near the stable man- 
ifold, after a transient time a steady state is reached where 
S-particles are accumulated along a fattened-up copy of the 
unstable manifold, as shown in Fig. ^d) for massless point 
particles, in Fig.^e) for bubbles, and in Fig.^f) for aerosols. 
In the computation, particles A and B are set to have the same 
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mass ratio and inertia parameters. To compute the effective 
fractal dimension D^f j of the unstable manifold, we use the 
uncertainty algorithm ITtIi applied to the first order approx- 
imation of the inverse map. The effective dimension turns 
out to be constant over many orders of magnitude of varia- 
tions in e and it is approximately the same for both nonin- 
ertial and slightly inertial bubble particles (De/f = 1.73 for 
e > 10^^^), while it is somewhat smaller for slightly inertial 
aerosol particles (Deff = 1.68 for e > 10"^'^), as shown in 
Fig.|3a). Strong evidence of the scaling law ^ is presented in 
Fig.|2jb) for two different values of the time lag r, where the 
scaling exponent is consistent with -De// = 1-73 for nonin- 
ertial and bubble particles, and with D^ff — 1.68 for aerosol 
particles. We see that even though the area Ab changes with 
the inertial properties of the particles, the scaling of Ab re- 
mains essentially the same for bubbles, as expected from our 
Cantor-set model. 

It is instructive to compare this result with the reaction-free 
dynamics. The dimension of nonhyerbolic invariant sets can 
be argued to be integer by mean of a zoom-in technique, where 
a fast numerical convergence is achieved by focusing on the 
densest parts of the fractal fTlll . The reaction, however, has a 
global character, as it takes place along the unstable manifold 
around the entire mixing region. This makes the convergence 
of the relevant effective dimension extremely slow, and that is 
why the effective dimension is apparently constant. 

In summary, we have shown that the dynamical system ap- 



proach to the reactive dynamics in imperfectly mixed flows 
also applies to realistic situations where nonhyperbolic and 
inertial effects are relevant. The rate equations of reactive 
processes are primarily governed by finite-time dynamics and 
as such change smoothly in the noninertial — > inertial transi- 
tion, which is in sharp contrast with the metamorphosis under- 
gone by the long-term and asymptotic dynamics of reaction- 
free particles. We have focused on autocatalytic reactions, 
but our results are expected to hold whenever the reaction 
front mimics the underlying unstable manifold and the activ- 
ity takes place along the boundary of a fattened-up fractal. 
Examples of this kind of process include infective reactions 
(e.g., combustion LI 8.1 ) and collisional reactions in general 
(e.g., A + B ^ 2C, where an unlimited amount of mate- 
rial A is present and material B is continuously injected in the 
vicinity of the stable manifold |3]). Finally, we observe that 
our analysis does not rely on the existence of a well-defined 
fractal set in the advection dynamics. The results remain valid 
as long as effective values for the fractal dimension and escape 
rate can be properly defined and are approximately constant 
over the relevant interval of observation. This is important for 
environmental processes, whose underlying dynamics is only 
partially understood fl3ll . 
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FIG. 1: For A = 4, (a) KAM island (light grey), and stable (grey) 
and unstable (black) manifolds for a = oo. (b) Fixed point attractor 
(black dot), basin of attraction (grey), and unstable manifold (black), 
for a = 1.05 and a = 1. (c) Stable (grey) and unstable (black) 
manifolds for a = 0.95 and a = 1, outside the region covered by the 
almost trapped orbits (light grey). Particles are launched with initial 
velocity matching the fluid velocity {5q — 0). (d-f) Corresponding 
area covered by B-particles in the "open" part of the flow right before 
the reaction, for r = 5 and a = 5.10"^. 



FIG. 2; For A = 4, (a) Effective dimension of the unstable manifold 
as computed from the uncertainty method, where /(e) is the fraction 
of £-uncertain points in the line a; = 0, < y < 0.1 of the time 
reversed dynamics, (b) Scaling of the relative area Ab covered by 
B-particles [in the region shown in Figs. Qd)-(f), right before the 
reaction] as a function of the reaction range a for two choices of the 
time lag r. In both plots, stars correspond to noninertial particles 
(a = oo), circles to bubble particles with a — 1.05 and a — 1, and 
plus signs to aerosol particles with a — 0.95 and a — 1. The aerosol 
data in (a) are shifted vertically downward for clarity. 
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